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<N : Abstract 

Although one-loop calculations provide a realistic description of bulk and single-particle nuclear 
properties, it is necessary to examine loop corrections to develop a systematic finite-density power- 
counting scheme for the nuclear many-body problem when loops are included. Moreover, it is 
■ imperative to study exchange and correlation corrections systematically to make reliable predictions 

for other nuclear observables. One must also verify that the natural sizes of the one- loop parameters 
are not destroyed by explicit inclusion of many-body corrections. The loop expansion is applied 
to a chiral effective hadronic lagrangian; with the techniques of Infrared Regularization, it is 
possible to separate out the short-range contributions and to write them as local products of fields 
that are already present in our lagrangian. (The appropriate field variables must be re-defined 
at each order in loops.) The corresponding parameters implicitly include short-range effects to 
all orders in the interaction, so these effects need not be calculated explicitly. The remaining 
(long-range) contributions that must be calculated are nonlocal and resemble those in conventional 
nuclear-structure calculations. Calculations at the two-loop level are carried out to illustrate these 
techniques at finite densities and to verify that the coupling parameters remain natural when fitted 
to the empirical properties of equilibrium nuclear matter. 
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I. INTRODUCTION 



Quantum Hadrodynamics (QHD) is a low-energy effective theory of the strong interaction 
that describes the strong nuclear force by the exchange of mesons between nucleons (the 
observed degrees of freedom at this energy scale). Initially proposed by Walecka [1], it 
has evolved over the years into a framework based on effective field theory (EFT) and 
density functional theory (DFT) [2-8]. EFT embodies basic principles that are common to 
many areas of physics, such as the separation of length scales in the description of natural 
phenomena. In EFT, the long-range (nonlocal) dynamics is included explicitly, while the 
short-range (local) dynamics is parametrized generically; all of the dynamics is constrained 
by the symmetries of the interaction. DFT tells us that the nuclear many-body system can be 
described by a universal energy functional that depends on nuclear densities and four- vector 
currents. With knowledge of the energy functional, one can calculate any observable for the 
(zero-temperature) many-body system. Moreover, a simplified treatment of the functional 
based on quasi-particle orbitals still provides an exact description of the bulk properties 
and some single-particle observables. Thus knowledge of the many-particle wave function is 
not needed to calculate this subset of observables [6]. The energy functional is constructed 
as an expansion in small parameters (the mean meson fields, which are in fact Kohn-Sham 
potentials, divided by a heavy mass, which could be the nucleon mass or the chiral symmetry 
breaking scale) and includes all possible terms consistent with the underlying symmetries of 
the system. 

In the current formalism, field redefinitions are employed to place the complexity of the 
problem in the meson field self-interactions. Each term is characterized by an unknown 
coefficient which, once all the dimensional and combinatorial factors have been removed, is 
a dimensionless constant of order unity, an assumption known as naturalness [9, 10]. The 
natural separation of length scales of the system is embodied in this theory: the long-range 
dynamics is included explicitly and the short-range physics is contained in the parametriza- 
tion. While this theory in principle contains all possible terms consistent with the underlying 
symmetries of QCD, in practice this is a perturbative expansion for the energy functional 
that can be truncated at a manageable level. Once it has been truncated, the now finite 
number of coefficients are fixed by experimental data; this theory can then be used for 
predictive purposes [11-14]. 

QHD is a strong-coupling theory. Unlike QCD, QHD is not asymptotically free, and 
currently no lattice version exists. Moreover, unlike QED, the QHD couplings are large, 
and there is no obvious asymptotic expansion to use to obtain results and refine them 
systematically. It is thus unknown whether QHD permits any expansion for systematic 
computation and refinement of theoretical results. One possibility is the loop expansion, 
which was partially explored in [15, 16]. In that work, explicit calculation of the short-range 
dynamics in terms of nucleons and heavy mesons produced enormous contributions to the 
energy, which rendered the loop expansion useless. This problem was only partially avoided 
by the inclusion of vertex corrections in the loops. In contrast, with the ideas of EFT, we can 
now provide a straightforward, physically motivated discussion of the short-range dynamics. 

In the present work, the loop expansion for QHD is constructed in the usual manner 
[15, 17-19]. First we define the action and the exact ground-state generating functional 
through a path integral. This generating functional contains all possible diagrams. The path 
integral is used to define the effective action, which is then expanded in powers of h (where 
h acts as a bookkeeping parameter and is not necessarily small [20-22]). An equivalent way 



2 



to state this is that the effective action is expanded around its classical value by grouping 
terms according to the number of quantum loops in their corresponding diagrams. Then 
we perform the functional derivatives and acquire the loop integrals. For the purposes of 
this work, we are interested in the expansion only up to the two-loop level. Consideration 
of the effects of contributions from the three-loop level and higher will be considered in 
future investigations [23, 24]. All of the integrals that represent tadpole and disconnected 
diagrams cancel out in the two-loop effective action, and we are left with only the fully 
connected diagrams. In the nuclear matter limit, the effective action is proportional to the 
energy density. For the cases considered here, there are three integrals of interest. These 
three integrals each have two factors of the nucleon propagator and one meson propagator 
(either scalar, vector, or pion). 

Why a loop expansion? The loop expansion is a simple and well-developed expansion 
scheme in powers of H that is derived from the path integral. The mean meson fields are 
included non-perturbatively and the correlations are included perturbatively. Therefore, one 
can analyze the many-body effects order by order. Indeed, previous work has shown the 
importance of "Hartree dominance" [25] : the mean-field terms dominate the nuclear energy, 
and exchange and correlation effects do not significantly modify the energy or nucleon self- 
energies, at least for states in the Fermi sea. We stress that we are not certain that the loop 
expansion is practical; the answer to this question is left for future consideration. However, 
the loop expansion has the advantage that it is fairly easy to separate the short-range and 
long-range dynamics and to analyze their structures. 

The nucleon propagator can be separated into two components, known as the Feynman 
and Density parts [19]. This is accomplished by taking into account the proper pole structure 
of the propagator. The Feynman part describes the propagation of a baryon or antibaryon; 
the Density part involves only on-shell propagation in the Fermi sea and incorporates the 
exclusion principle. 

As a result, the two- loop integrals can each be separated into three distinct parts, which 
we refer to as exchange, Lamb-shift, and vacuum-fluctuation contributions [15]. The ex- 
change term has two factors of the Density portion of the nucleon propagator. Thus both 
momentum integrals are entirely within the Fermi surface and the exchange term is finite. 
This represents a contribution from long-range (nonlocal) physics and must be calculated 
explicitly. The Lamb-shift term contains both Feynman and Density parts of the nucleon 
propagator. This term is so named because it is analogous to the Lamb shift in atomic 
physics, where a particle in an occupied state interacts with a virtual particle in an unoccu- 
pied state that shifts its spectrum. This contribution is short-range, and we will show that 
it can be expressed as a sum of terms that are already present in the QHD EFT lagrangian. 
As this is an effective theory, the coefficients of these terms are determined by matching 
to empirical data; thus, these short-range terms are just absorbed into local terms already 
present in the lagrangian and should not be calculated explicitly. The vacuum fluctuation 
term contains two factors of the Feynman propagator and both momentum integrals extend 
outside the Fermi sphere. This involves the excitation of NN pairs and is therefore also 
short-range physics. We will show that it can also be expressed as a sum of terms which 
exist in the EFT lagrangian. As a result, it can be removed in the same manner as the Lamb 
shift. In addition, if nonlinear meson self- interactions are included, a number of pure meson 
loops arise. These terms, however, can be expressed as a power series in the meson fields 
with undetermined coefficients. As before, these terms are just absorbed and do not need 
to be calculated explicitly. The result is that, for the cases considered here, only three finite 
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integrals representing the two-loop contributions from the scalar and vector meson and the 
pion need to be calculated. 

It is interesting to note that the procedure described above is similar to Infrared Regular- 
ization. In Infrared Regularization, the one-loop self-energy contribution can be separated 
into soft and hard parts [26-30] (they are sometimes referred to as the infrared singular and 
regular parts [26, 27]). The hard (or infrared regular) parts arise from large momentum, or 
high-energy, contributions and are expressible as a power series of terms already contained 
in the underlying lagrangian; as in our case, they are just absorbed into the coefficients. The 
soft (or infrared singular) parts, in the notation of Ellis and Tang, contain both analytic and 
nonanalytic terms. The analytic portion results from high-energy dynamics and is therefore 
treated in the same manner as the hard contributions. The nonanalytic portion develops 
from low momentum, or low-energy physics, and is essentially long-range dynamics. This is 
the nonlocal contribution that must be calculated explicitly. The regularization procedure 
for the (closed) energy loops in this work is significantly simpler than the procedure for 
diagrams with external momenta, like scattering amplitudes or self-energies, because the 
momenta carried by the boson propagators are always spacelike. 

In the QHD lagrangian, a well-developed mean-field power-counting scheme has been 
devised [2, 3]. As noted, the energy functional is an expansion in small parameters (the 
meson fields divided by the nucleon mass). In addition, the Fermi wave number, which is 
related to the size of the derivatives, is also small when divided by the nucleon mass. If 
naturalness holds, then adding up the powers of these ratios (and some counting factors) 
yields an accurate estimate of the size of a given term. This paper will seek to investigate 
the relationship between the two-loop integrals and this underlying power counting scheme. 

The purpose of this work is to illustrate these techniques by working to the two-loop level 
in a loop expansion; the resulting two-loop integrals are separated into long-range and short- 
range physics. The short-range contributions are expressed in forms that already appear at 
the mean field level. Since the coefficients of these terms have yet to be determined, they are 
just redefined, thereby incorporating these new contributions into the mean field lagrangian. 
As a result, they are already present in the one- loop level QHD calculation. The long-range, 
nonlocal physics must be explicitly calculated. In this work, we fit the two-loop energy to 
the equilibrium point of nuclear matter. Then we compare sets constructed at the two- 
loop level with those of previous work developed at the mean field level and consider the 
naturalness of the parameters. We also examine how these new contributions fit into the 
power counting scheme. Some of this work has been discussed previously in an unpublished 
Ph.D. dissertation [31]. 

II. THEORY 

In this section, we follow the usual procedure for constructing the loop expansion [19] 
and in the following subsections examine the one- and two-loop contributions, as in [15]. 

A. Loop Expansion — Background 

Consider the following nonrenormalizable effective lagrangian, which extends the Walecka 
model to include the nonrenormalizable nN coupling and nonlinear scalar and vector meson 
interactions. This lagrangian can also be obtained from the chirally invariant lagrangian 
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in [2] by retaining only the lowest-order terms in the pion fields — sufficient for the two-loop 
calculations in this work — and a subset of the meson nonlinearities 1 : 



C = -i> 



7m {d» - igvVn) 



i^7~lvtlsdvtE. + (M — g.s4>) 

Jit 



+ c NL + sc , 



(1) 

where V^„ = d^V v ~d u V^ and 7r = jn a -r a . Here ip are the fermion fields and 0, V^, and n a are 
the meson fields (isoscalar-scalar, isoscalar-vector, and isovector-pseudoscalar, respectively) 
and the heavy meson fields are also chiral scalars. For the purposes of this work, we will not 
include effects from the isovector-vector channel or the electromagnetic field, but to do so 
is straightforward [2, 32]. The numerically small tensor coupling between the omega meson 
and the nucleon is not included, and nonlinearities in the meson sector of the lagrangian 
(Cnl) will be considered later at the mean field level; the study of their effects at two-loop 
order is left for future work [23]. S£ contains all of the counterterms. Note that in this 
work, the conventions of [8] are used. 
The action is defined as 

Sl0,tg = J d 4 xC(x) , (2) 

and the exact ground-state to ground-state generating functional is defined through the 
following path integral [15] 

Z[j,J^] = ex P {tW{j,J,]/h} 

= M- 1 J DfflDWDWDiVjDM 



x exp 



h 



d 4 x [C{x)+j{x)(t>{x) + J lx {x)V ll {x)\ 



where 



D{^)D^)D{<j>)D{V ll )D{'K a )^\^ / d 4 xC(x) 



(3) 



(4) 



Here M is the normalization factor (in effect, the vacuum subtraction), and j(x) and J^x) 
are the external sources corresponding to the meson fields and V^, respectively. The 
connected generating functional W\j, J M ] contains all of the connected diagrams, and Z[j, J M ] 
contains all possible diagrams without "vacuum bubbles" . 

The classical values of the meson fields are determined by extremizing the action: 



5S 

5<j){x) 



= (d 2 - m 2 s ) 0o = -j(x) 



4>=<f>o 



ss 



SV^x) 



m 



v 



V° = -Ux) . 



(5) 
(6) 



1 We are not making a chiral expansion in powers of the pion mass. We have simply included a pion mass 
for kinematical purposes in Eq. (1). 



5 



Next, we replace the fields by quantum fluctuations around their classical fields: ip(x) — > 
^ 1/2 V>(x), ip(x) -> h l / 2 %l)(x), (f)(x) -> 0o(a?) + h^ 2 <r(x), V^(x) -> V^z) + h^fj^x), and 
7i"a(^) - ► ^ 1 ^ 2 fi a (^)- Notice that the fermion fields have no classical limit (and the pion 
has no mean field if one assumes no pion condensate). Observe that a factor of h 1 ^ 2 is 
associated with each quantum fluctuation; the numerical value of h is immaterial, as it is 
just a bookkeeping parameter [20-22]. Next, a number of extra sources are included that are 
set equal to zero at the end (u, U^, ( a , £, and £ corresponding to the quantum fluctuations 
a, fjfj,, Q a , tp, and ip, respectively). The path integral is then rewritten in terms of functional 
derivatives with respect to these new sources. The final result is [15] 



Z\j,J li ]=M'- 1 exp{- I d 4 x [C (x)+j(x)Mx) + UxK(x)} 



x exp {tr In [G^Gjj 1 ] } 



45 



exp \ ih 1 ^ 2 I d 4 x 
-i6 



iS 



6t(x) 



-iS 



-W +9s 



5u(x) 



-iS 



x exp < —i 



x 




[Ka{x)\ 2 

d A xd A yi{x)G H {x-y)i{y) 
exp |^ J j d A xd A y [u(x)A° s (x - y)u(y) 
+ U,(x)V%(x - y)U v (y) + Ca(x)Af(x - y)( b (y)] 



6t(x) 



(7) 



-J sources = 



where A/ 7 is equal to the portion of Z involving the variational derivatives, but with free 
propagators. £ represents the lagrangian at the mean field level. The fermion propagators 
in momentum space are 



G° F (k) 



G H (k) = - 



ift + M — ie ' 
-1 



i# - + (M - g S (f>o) ' 

and the free meson propagators in momentum space are 

1 



(8) 
(9) 



A°(fc) = 

KM = 

K\k) = 



k 2 + m 2 s — ie ' 
1 

k 2 + my — ie 



r k^k v 



m 



V 



1 



k 2 + ml — ie 



>ab ■ 



(10) 

(11) 
(12) 



It is assumed that all of the divergent integrals are regularized in some fashion that 
preserves the symmetries of the theory, for example, dimensional regularization. The lon- 
gitudinal term in the vector meson propagator vanishes in the following analysis, as the 
vector meson couples to the conserved baryon current [19]. As usual, the infinitesimal e is 
introduced to generate the proper pole structure. 



6 



The expectation values of the meson fields in the presence of external sources are 



= (O + |0(x)|O-) = _ ih SlnZ[j, J,} _ 5W[j, J,} 



(0+|0-> 
(0+|Vi(x)|0-> 



Sj(x) 



Sj(x) 



_. h 5\nZ\j,J,} _ 5W\j,J^} 



(0+|0-) SJ^x) 5Ux) 

In the limit of vanishing sources, these expectation values become 

lim e = = constant , 
lim Vf, — V u — constant . 

We now define the effective action T by a functional Legendre transformation: 

r[0 e , v;\ = w[j, j,] - J d 4 x [j(x)<j> e (x) + U^Kix)] . 

In uniform nuclear matter, this effective action is related to the energy density S by 



lim r[0 e ,y;] = - / d^s^v,} 



(13) 
(14) 

(15) 
(16) 

(17) 

(18) 



B. One-Loop Calculation 

In this section, we consider the lowest-order terms in the loop expansion, the one-loop 
contributions; these are the terms that correspond to 0{h°) in the expansion. The generating 
functional at one-loop order is [15] 

zW[j,J^\ ee exp{tW^[j,J,}/h} 

= expj^ J d 4 x [C (x)+j(x) ( j )0 + .Ux)V^ 

xexpjtr ln[G°G#]} , (19) 

where we have assumed uniform classical fields. Therefore, the connected generating func- 
tional is (for now we have dropped the nonlinear meson self-interactions) 



d 4 x l--m 2 ^ 2 ^t/Ot/O 



f d 4 k 
—in I - — — 

J (2ir) 4 



tr In 



ift + M — ie 



(20) 



Here "tr" denotes the summation over both spin and isospin. Following from the previous 
section, Eqs. (13) and (14) become (f)p(x) = 0o and V^ 4 \x) = V®, respectively. For a spa- 
tially uniform system, the classical vector field is V® = iVoS^. To calculate the momentum 
integral, we use the relation [15] 

d 4 k 



-in 



ih 



(2^) 4 
d 4 k 



tr In 



1 + 



gV74Vb - g,s4>o 



tr[ l4 h(G H (k)-G° F (k))] , 



(21) 
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where G F comes from the vacuum subtraction A/", which determines the proper bound- 
ary conditions. Then, using dimensional regularization to eliminate the vector field in the 
momentum integrals, the effective action becomes 

r (1) [0 e , V e °] = J d 4 x \^m 2 v V 2 - \m\<g - hg v V oPB 

+ih J-0ji k * tT faGM*)] - ihj -J^ htr [^G F (k)} } , (22) 

where E*(k) = (k 2 + M* 2 ) 1/2 , the effective nucleon mass M* = M - gs<fio, and 7 is the 
spin-isospin degeneracy. Here we have separated the nucleon propagator into two parts [19] 
(known as the Feynman and Density contributions), as shown below: 



G*(k) ' ' 



= (i]6-M* 



?7T — * 

-5[k 4 -E*(k)]6(k F -\k\) 



_k 2 + M* 2 -ie E*(k) 
= G* F (k) + G* D (k) ■ (23) 

This is done by taking the proper pole structure into account. The Feynman part describes 
the propagation of baryons and antibaryons; the Density part involves only on-shell propa- 
gation in the Fermi sea and corrects the propagation of positive-energy baryons for the Pauli 
exclusion principle. 

To remove the divergences in the momentum integral in Eq. (22), we expand G* F (k) as a 
polynomial in gs4>o (the Furry expansion [15]) 



G* F (k) = J2(M*- Mf [G° F (k)} n+1 + (AT - M) m+1 [G° F (k)] m+1 G* F (k) , (24) 

?1 = 



which is valid for any m > 0. If we take m = 4 and insert into Eq. (22), we notice that the 
first four terms are divergent. The counterterms introduced to absorb these divergences are 

4 

5CM = hY,<*n4>o, (25) 

n = 1 

where aij are pure numbers. The first term in the Furry expansion cancels the last term in 
Eq. (22); the next four terms cancel with the counterterms in Eq. (25). Thus, the one- loop 
effective action is just 

[0 e , K°] = / d'x I \m 2 v V 2 - l -m 2 4>l - g v V p B 

~J2^ J d " kE ^ k ) 9 ^ - \ k \) - A£W(M*)j , (26) 
where the factor h has been omitted and 

A£ VF {M*) = -i{M*-Mf j 0- 4 k 4 tr{ l4 [G F (k)fG F (k)} . (27) 
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As before, the mean fields are determined by extremizing the effective action. At the 
one-loop level, the energy density is [using Eq. (18)] 



£U[M\p B \ = g v V oPB - l -m 2 v V* + ^|(M - M*f 

+ (2^ j d3kE *^ kp ~ ^ + A£ vf(M*) , (28) 



which is the relativistic Hartree approximation in the original Walecka model [1]. The final 
term in the energy density is written as [15] 

4!(n-5)!(^o) w 1 , 9Q . 

Note that this term has no explicit density dependence. 

As discussed in [3, 33], the general form of each term in A£ VF shows that this vacuum 
contribution is unnaturally large. If we accept the naturalness assumption, the conclu- 
sion is that the vacuum contribution is not well described by ASvf at the one-baryon-loop 
level. Baryons are incorrect degrees of freedom for computing short-range loops. Further- 
more, to make these terms natural, it is certain that there must be large cancellations from 
higher-order loops, which therefore must be calculated. Yet, we know that in principle the 
polynomial terms in should appear in an effective lagrangian, since they satisfy all the 
symmetry requirements; moreover, there will always be short-range contributions to these 
terms that we cannot calculate. So the proposal is that we need not work hard to get the 
vacuum contributions — instead, we can adjust the unknown natural coefficients to include 
them! 



C. Two-Loop Calculation 



In this section, we present the corrections to the theory arising from the two-loop contri- 
butions. We define the connected generating functional at the two-loop level as 



W (2) = W (i) + w 



2 ■ 



(30) 



Keeping only the terms of 0(h) in the loop expansion (which is essentially an expansion in 
the coupling constants gs, fiV, and g^), the connected generating functional is 



W 2 



2 




d xd y\g s 5 afS 5 a/ p 



~ 9v (l^)ap(lu)a'P' 
7^75*9, 



-iS 



SU^x) 



-iS 
Su(x) 
-iS 



-iS 
5u(y) 
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i5 



Ka{x) 



2 



-iS 



n 
2 



a'P' 



9 



X 



" iS ' 



' -iS ' 

a M*). 



' iS ' 

Mv). 



' -iS ' 
My). 



expj-i^ J d 4 xd 4 y£(x)G H (x - y)£(y) 
d 4 xd 4 y [u(x)A° s (x-y)u(y) 



x 



x exp <; - 



+ U a (x)V° ap (x - y)U p (y) + ( c (x)A c *(x - y)Uv)] 
-VEV , 



sources = 



(31) 



where VEV is the vacuum subtraction, which is just equivalent to the rest of Wi with free 
propagators. After working out the variational derivatives, the full expression becomes 



Wo 



h 2 



d x 



d A k d 4 q 
(2tt) 4 (2tt) 4 



x ( gl { A° s (k - q) tr [G H (k)G H (q)\ - A° (0) tr [G H (k)\ tr [G H (q)\ } 

9v {VU k - 9) tr [l,G H (k) lu G H (q)] - V%(0) tr [^G H (k)\ tr [ lu G H (q)]} 



9 A a ah 



-VEV 



Af(k - q) tr 



Al B jG H {k)(f6- 4)l^G H (q) 



(32) 



W 2 contains all the connected two-loop diagrams. To isolate and remove the tadpole 
diagrams, we must consider the effective action. The expectation values of the meson fields 
at the two-loop level are [15] 



<t>f\x) 



6w<u 

5j(x) 
8WV> 
5J^(x) 



(x) -hg s d 4 y 



V°Jx)-ihg v / d 4 y 



d 4 k 

d 4 k 



(2k) 



A° s (x-y)tr[G H (k)} , 



4 V%(x-y)tr[ ltM G H (k)} , (33) 



where = S + W 1: and there are no pion tadpoles. Thus, the effective meson fields 

contain tadpole contributions. 

The effective action at the two-loop level is 

+ J d 4 x {j(x) [Mx) - Ux)\ + W [v°(x) - v;{x)] } 

+ 0(h 3 ) . (34) 

Next, we change variables to O = 4>e + 4>i an d — + Vp, where <j>i and cancel 
quantum corrections contained in <p e and respectively. Then, we expand this expression 
as a power series in <f>i and about <p e and V*. Using Eqs. (5) and (6), we get 
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ft' 2 , 

+ y9§ 



d x 



d A k d A q 
(2tt) 4 (2tt) 4 

d 4 k d A q 
(2tt) 4 (2tt) 4 



A°(0)tr[G H (A;)]tr 

2^(0) tr [ 7# ,Gh(*)] tr[ 7 ^( g )] , (35) 



where the last two terms cancel with terms in W2 (the tadpole diagrams drop out). Now 
we use dimensional regularization to eliminate the dependence on V®. Lastly, we write the 
energy density as £^ = £^ + £2, where [15] 



So 




d A k d A q 
(2tt) 4 (2tt) 4 



9 s ao 



A° s (k - q) {tr [G*(k)G*(q)\ - tr [G F (k)G° F (q)} } 



9 f V%(k - q) {tr [ lfl G*(k) lu G*(q)] - tr [ lfl G F (k) lu G F (q)} } 
~ ^-Afik - q) {tr p- ^G*(m- 4)^G\q)] 



tr 



(ft- 07fcyG° (*)(*- d_h^G F (q)]} 



(36) 



and the factors of H have been suppressed. The terms involving the free nucleon propagators 
come from the VEV subtraction. The corresponding two-loop diagrams are shown in Fig. 1. 

Since each G* can be separated into Feynman and Density parts [see Eq.(23)], we can 
rewrite the first term in £ 2 as the following sum [15] 



£ 



(2) 



c(2) , M2) p {2) 

^ -EX ' ^rh—T.G T 



-LS 



VF 1 



(37) 






FIG. 1: Two- loop diagrams. The double lines represent baryon propagators [Eq. (23)]. The dashed, 
wiggly, and dotted lines represent scalar, vector, and pion propagators, respectively [Eqs. (10) to 
(12)]- 



where 



r(2) 
'4>-EX 

c(2) 
C 4>-LS 



,(2) 

'4>—VF 



9l 



-9 2 s 
9l 






d A k d A q 
(2tt) 4 (2tt) 4 
d A k d A q 
(2tt) 4 (2vr) 4 
d A k d A q 
(2tt) 4 (2tt) 4 
x{tr[G F (k)G F (q)]-tr [G° F (k)G° F (q)] } 



A^-^tr^)^)] , 

A^-^tr^)^)] , 
A° (k - q) 



(38) 
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It turns out this is also true for the vector and pion terms, so that we may write the energy 
density at the two-loop level as 

SV = £U+SV+S%+£V. (39) 

(2) 

£ EX is called the "exchange term". It has two factors of G* D , which restricts the double 
integral to the inside of the Fermi sphere. It is finite and explicitly density-dependent; 
therefore, we regard it as a purely many-body effect. It can be calculated straightforwardly. 
This corresponds to the exchange of identical fermions in occupied states. 

(2) 

S L g is analogous to the Lamb shift in atomic physics, since it involves particles in occu- 
pied states whose energies are shifted by interactions with the fluctuating meson fields at 
finite density. These fluctuations modify the short-range structure of the baryon due to the 
existence of the background meson fields. 

Syp is a true vacuum fluctuation, since it involves both meson and baryon virtual exci- 
tations. This involves the excitations of NN pairs, which is also short-range physics. 

(2) 

The £ FX contributions contain long-range dynamics [they are characterized by the length 
scale of 0(kp 1 ) ~ 0.8 fm] and are nonlocal (they involve logarithmic functions of momenta); 
thus, we must calculate them explicitly. The contributions £ E g and £y F correspond to short- 
range physics [they are characterized by length scales of O^M" 1 ) ps 0.2 fm]. So we expect 
that the latter two can be absorbed into the coefficients of the local terms in the lagrangian. 

Normally, certain counterterms are introduced at the two-loop level to deal with the 
divergences in £^s and £y F , and the remaining finite contributions, which depend on the 
renormalization conditions, can be calculated numerically [15]. However, for a nonrenor- 
malizable effective lagrangian, an infinite number of counterterms are required in principle. 

(2) (2) 

Here we will argue that both £ L $ and £ VF can be written in forms that are already present 
in our lagrangian (before truncation). Therefore, by adjusting the coefficients of these terms, 
these contributions are completely absorbed. 



D. Short-range Dynamics 

Consider the Lamb-shift contribution from the scalar meson 

d 4 k d 4 q 



,(2) 
■<h-LS 



-9l 




(2tt) 4 (2tt) 



A° s (k-q)tr[Cr F (k)Gh(q)] 




d 4 k 



tr 



where the self-energy is 



9l 



d 4 q 



A° s (k - q)G* F (q) . 



(40) 



(41) 



Note that this expression contains no explicit density dependence. We substitute the Furry 
expansion [see Eq. (24)] into the self-energy: 



2 00 r j4 

= -f £ ( M * - M T j (^4 A °s( k - 9) [<%)] 



n+1 



(42) 
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but now we have let the sum go to infinity and not cut it off as before. This can be rewritten 
as a Taylor series 



n = 



(m* - My 



(43) 



M*=M 



We then expand Y?p{k) as a Taylor series around ift = M to separate the momentum 
dependence, or 



E^(Jfe) = Y#(iJ6 = M) + 



d(i#) 



1 d 2 E^(A;) 



■;^=m 



2 



(i#-M) 2 + ... 



(44) 



The self-energy now becomes 

F V ^ m ! n ! d(ij() m dM* n 



m=0 n=0 

x (i# - M) m (M* - M) n 



i/i=M,M*=M 



oo oo 



E E w - (m* - m)» , 



(45) 



m = n = 



where a^ n are just pure numbers, some of which are finite and others infinite (these corre- 
spond to divergent diagrams). This consequence of Lorentz covariance allows one to rewrite 
the Lamb-shift contribution to the energy density as 



oo oo 



,(2) 



-LS 



= ^E E 



m = n = 

oo oo 



d A k 



= 2 E E°-(M*-Mr / 

(k F - |jfe|) 



tr [G* D (k) {if- M) m (M* - M) r 



(2^) 4 £*(fc) 



m = n = 

x5[fc 4 -£*(£)] (A*- |Jfe|)tr[(i#- M*)(i#-M) m ] . 



In general, one can write 

(if- M) m = f m (k 2 , M) if + g m {k\ M) , 
where f m and # m are polynomials in k 2 and M. This lets one reduce the integral to 

d A k 



(46) 



(47) 



-2 7 



x 



(2tt)< 



[A; 2 / m (A; 2 ,M) + M*£ m (A; 2 ,M)] 



27T 



E*(fc) 



S[k 4 -E*(k)}9 (k F - \k 



= - F m (M*, M)—^-— 
2 v ' ; (2tt) 3 

= ±F m (M*,M) Ps . 



d 3 k 



M* 



E*(k) 



kp - \k\ 



(48) 



13 



Here ps is the scalar baryon density and the on-shell condition has been used. Substituting 
this result into Eq. (46), we get 



,(2) 
'6—LS 



PsJ2 E a tnFm(M* , M) (M* - MY 

m=0 n=0 

oo 

Ps £ ai(gs<P o r ■ 



(49) 



rn = 



Before the field redefinitions were conducted on the effective lagrangian to put it in canonical 
form, terms like ps{gs4>o) n — (V^Hfi's^o)™ appeared in the theory Since the coefficients of 
these terms will be eliminated by field redefinitions, the two-loop contributions to these terms 
can just be absorbed, and there is no need to calculate them. These arguments rely only on 
the Lorentz structure of the self-energy and the on-shell condition imposed by the Density 
propagator, so similar conclusions follow for the contributions from the vector mesons and 
pions. 

We treat the term £j?l VF in a similar manner: 



,(2) 

'S-VF 



9l 



d 4 k d A q 
(2vr) 4 (2tt) 4 



A° s (k - q) {tr [G* F (k)G F (q)\ - tr [G F (k)G F (q)] } . 



(50) 



Note that this contribution has no explicit density dependence. Using Eq. (24), we can write 

d 4 k d A q 



,(2) 



9 2 s 



-VF 




(2tt) 4 (2tt) 4 



A° s (k - q) 



oo oo 

E E 

m= n = 



(M* - M) 



rn+n 



X 



1 


m+l 


1 


k 2 + M 2 _ 




_q 2 + M 2 _ 



n+l 



tr [{i}t-M) m+1 {i4-M) n+1 ] , (51) 



where m + n^0. (The m + n = term cancels the vacuum term.) We can use Eq. (47) to 
solve the trace, or 



£ 



(2) 



-VF 



-ah 




d 4 k d 4 q 



x 



(2tt) 4 (2tt) 4 

oo oo 

£(M*-M) m+n 



A° s (k - q) 



1 


m+l 


1 


n+l 


_k 2 + M 2 _ 




_q 2 + M 2 _ 





m = n = 

x [g m+1 (k 2 ,M)g n+1 (q 2 ,M) - k ■ qf m+1 (k 2 , M)f n+1 (q 2 , M)} 

,2 _ j\/r2 „„j „2 _ T\^2 



(52) 



Now one can expand the baryon denominators around k = M and q = M and perform 
a Wick rotation to Euclidean space. What remains are some integrals over a complicated 
polynomial in k 2 , q 2 , k-q, M*, and M. However complicated, these four- dimensional integrals 
can be done, and we are left with 



(53) 



m = 1 



where are constants that depend on M. These terms can also be absorbed into preexisting 
terms in the lagrangian, and hence there is no need to calculate them. Fortunately, only 
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a few coefficients are required for an accurate description of bulk nuclear properties [34]. 
Note that the forms in Eqs. (49) and (53) are consistent with the explicit results for these 
integrals given in [15]. 



E. Long-range Physics 



In this section, we consider the two-loop exchange contributions to the energy density 
in the presence of background mean fields. These terms are nonanalytic functions of the 
Fermi momentum and correspond to nonlocal contributions to the energy. The scalar meson 
contribution to the exchange term becomes 

d A k d A q 



?(2) 
'<fi—EX 



9l 




A»(k-q)tr[GUk)G* D (q)] 



9 2 s 




(2tt) 4 (2tt) 4 
d A k d A q 9{k F -\k\)9{k F 



151) 



Tgg 
32tt 4 



x 



(2tt) 3 (2tt) 3 2E*(k) 2E*(q) 
x5(k 4 -E*(k))S(q 4 -E*(q)) 
x ti[(ifl-M*)(ift-M*)] 
kF \k\ 2 d\k\ f kF \q\ 2 d\q\ 



(k — q) 2 + m| 



E*(k)E*(q) - |A;||g1cos^ + M 



*2 



(54) 



2E*{k)E*{q) - 2\k\\q\ cos6 - 2M* 2 + m% 

where we have used Eq. (23) and integrated out most of the angular dependence. The vector 
meson contribution is 

d A k d 4 q 



(2) 

V-EX 



9v 




9v 




d A k d A q 6{k F - \k\) 9(k F - \q\) 
(2tt) 3 (2tt) 3 2E*(k) 2E*(q) 



x5(h-E*(k))5(q 4 -E*(q)) 



(/.• q)' 2 ■ in' 2 



'v 



xtr [ llM {i^-M*) lv {i4-M*)\ 



19v_ 
16tt 4 



k\ 2 d\k\ 



q?d\q\ 



x 



, d(cosO) 
E*(k) J E*(q) J_ x 1 J 

E*(k)E*(q) - \k\\q\ cosfl - 2M* 2 
2E*(k)E*(q) - 2\k\\q\ cos6 - 2M* 2 + m 



v 



(55) 



where we integrated out most of the angular dependence. The contribution to the integral 
from the longitudinal portion of the vector propagator vanishes if one works out the trace. 
Finally, we consider the two-loop contribution arising from pion exchange. We remove most 
of the angular dependence and arrive at the result 

d A k d A q 



r(2) 
y ix—EX 




(2tt) 4 (2tt) 



Af (k - q) 
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x tr 



9l 




d 4 k d A q 9{kp — \k\) 9{kp — \q\) 



(2tt) 3 (2tt) 3 2E*{k) 2E*(q) 
x5(k 4 - E*(k)) 5 (g 4 - E*(q)) — — I 



+ Till 



x tr 



(^-^)75y(^-M*)(^-g0 75 



14- M*) 
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57 — 



at 



E*(k)E*(q) - \k\\q\ cosfl - M* 2 



d(cos9) 



x 



2E*(k)E*(q) - 2\k\\q\ cos6 - 2M* 2 + ml 



(56) 



III. DISCUSSION 

In this work, we performed the loop expansion to the two-loop level for a model QHD 
lagrangian. This expansion provided a simple scheme for separating the short- and long- 
range dynamics order by order and for analyzing their structures. The short-range, local 
physics (including all the divergences) was isolated and absorbed into the coefficients of the 
lagrangian. The long-range dynamics is nonlocal and has to be calculated explicitly. For 
the purposes of this work, we are interested in applying the loop expansion only to ordinary, 
symmetric nuclear matter. For that reason, we neglected both electromagnetic effects and p 
meson exchange, although they both appear in the mean field lagrangian [2, 3]. In addition, 
we exclude the tensor terms in the fermion sector and the nonlinear meson self-couplings. 
The exploration of the effects of these terms at the two-loop level is left for future work. The 
inclusion of the two-loop effects also calls into question whether the naturalness assumption 
[2, 9, 10] is preserved; we partially answer that question here, but a more complete proof 
that naturalness holds is also left for a subsequent investigation. 

Now we consider the numerical analysis of the three surviving two-loop integrals (the 
exchange terms). We use the parameter sets listed in Table I. For the pion term, we use 
g\ = 1.5876 and f w = 93 MeV [2]. In addition, both L2 and Wl sets lead to equilibrium 
at k-p = 1.3 fm _1 . The mean meson fields are determined by extremizing the meson field 
equations and are used as input to the exchange integrals. Then the full two-loop energy 
density is extremized with respect to the meson fields. The results for all three exchange 
terms are shown in Table II for both L2 and Wl parameter sets. The third set (MOA) listed 
in Table I was fit at the two-loop level to nuclear equilibrium (£/ p B — M = —16.10 MeV and 
k-p = 1.3 fm _1 ) by adjusting g$ and gy using a downhill simplex method to minimize a least- 
squares fit (with respective weights of 0.0015 and 0.002). The binding curves of all three sets 
are shown in Fig. 2. The curves for the L2 and Wl sets include the two-loop contributions 
and are compared to the saturation curve of the MOA set at the two-loop level. Fig. 3 shows 
results when the L2 and Wl sets are evaluated at the one-loop level, and the MOA set is 
evaluated at the two- loop level. First, we note that nuclear saturation can be reproduced at 
the two- loop level with a parameter set (MOA) that is natural (gs,v ~ 47r). Second, Fig. 2 
shows that while the two-loop contributions are not large, they are not negligible either. 

Next, we compare the magnitudes of these integrals to terms in the meson sector at 
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t 9 r^i 


wi roi 

VV 1 [ZJ 


i\/rn A 


m s /M 


0.55378 


0.60305 


0.54 


m v /M 


0.83387 


0.83280 


0.83280 




0.83321 


0.93797 


0.79361 


g v /4ir 


1.09814 


1.13652 


0.96811 



TABLE I: Parameter sets used in this work. 





L2 [3] 


Wl [2] 


MOA 


c(2) 

°4>-EX 


42.59 


46.65 


40.30 


c(2) 


-29.51 


-30.60 


-23.14 


c(2) 

C TT-EX 


12.39 


12.21 


12.44 



TABLE II: Size of two-loop integrals for the parameter sets in Table I. Values are in MeV. 
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FIG. 2: Comparison of the nuclear binding curves for the sets L2 and Wl (both with one-loop 
parameters but with the two- loop contributions included) and MOA (two-loop level). 



the mean field level. However, neither set L2 nor Wl includes nonlinearities in the meson 
fields. Therefore, we will compare the exchange integrals to terms from the sets Ql and Q2 
[2]. Here Ql and Q2 both include cubic and quartic scalar field self-couplings and Q2 also 
contains a quartic vector field self-coupling. These nonlinear terms are contained in 

r f K 3 gs<fi . ^A9 2 s4> 2 \ 2 i2 . 1 > 2/ t/t/ «n2 /, 7 x 

= - [uJa + 4, -Jp ) m s^ + 4i Co 9v (V,Vn (57) 
of Eq. (1). This does not imply that these nonlinearities were taken into account when the 
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MOA 

-- L2 




FIG. 3: Comparison of the nuclear binding curves for the sets L2 and Wl (both one-loop level) 
and MOA (two- loop level). 



two-loop integrals were calculated; it is simply instructive to compare the relative sizes of 
the mean field nonlinearities and the two-loop contributions. Note that the Ql and Q2 mean 
fields were obtained by extremizing the one- loop energy, including Cnl, while the M* in the 
exchange integrals was obtained by minimizing S^ 1 of Eq. (39). While the inclusion of the 
nonlinear terms in the meson sector at the two-loop level has yet to be performed, a cursory 
glance tells us that they will affect only the meson propagators in two-loop integrals. Thus 
the overall magnitude of these terms is not expected to change much when these nonlinear 
effects are included. 

Fig. 4 shows the results of this comparison. The crosses in Fig. 4 represent the expected 
magnitude per order of terms in the meson sector using the rules of naive dimensional 
analysis (NDA) [2] with the chiral symmetry breaking scale of A = 650 MeV. Thus, one 
observes that the two-loop integrals are roughly equivalent to v — 3 order in the power 
counting in the mean field lagrangian. While this is not large, they cannot be neglected 
in a description of nuclear matter saturation properties, particularly in view of the nearly 
complete cancellation of scalar and vector terms at order v — 2. We emphasize that the 
explicit computation of short-range, two-loop contributions leads to unnaturally large terms 
in the energy of nuclear matter [15]; in the EFT approach, these contributions are absorbed 
in the local counterterms and are either removed by field redefinitions or are determined by 
fitting to empirical nuclear data. 

We stress that the results shown here are for the crudest truncation in the underlying 
lagrangian. It is well known that an accurate description of nuclear saturation (and the 
equation of state) requires the inclusion of nonlinear meson self-interactions, which is the 
subject of future work [23]. Furthermore, the question of whether or not the loop expansion 
is valid for QHD cannot be answered at the two-loop level. While two-loop calculations are 
a necessary step in that direction, an investigation into higher-loop effects is required; this 
is also the subject of future work [24]. 
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FIG. 4: Comparison of the magnitudes of the mean field terms in the meson sector with the two- 
loop exchange integrals. The inverted triangles represent, from top to bottom, the scalar, vector, 
and pion two-loop integrals. The abscissa represents the order v in the power counting at the mean 
field level [2]. Absolute values are shown. 



A. Infrared Regularization 

The separation of the short-range and long-range dynamics, and the subsequent absorp- 
tion of the short-range physics into the lagrangian are analogous to Infrared Regularization 
in Chiral Perturbation Theory [26-30]. In the language of Ellis and Tang [26, 27], a loop 
integral can be separated into hard and soft components 

G — (G - RSG) + RSG , (58) 

where G is the loop integral, S projects out the soft part, and R renormalizes SG to re- 
move the ultraviolet divergences. To extract the soft part, do the following: take the loop 
momentum to be of order Q, make a Q/M expansion of the integrand, and interchange 
the order of integration and summation. Now consider a loop integral with momentum q, 
specifically the q^ part of the integral. Closing the contour by a semicircle at infinity, we get 
the sum of three contributions: the semicircle, the soft poles, and the hard poles. The soft 
and hard poles are of order Q and M respectively. (It has been assumed that the hard and 
soft poles can be separated, as is the case in theories with Goldstone bosons and massive 
baryons.) The semicircle may produce divergences, but these can be removed by the usual 
renormalization. The soft poles cannot be expanded in Q/M, such as those in the pion 
propagator. However, a Q/M expansion can be made around the hard poles since the loop 
momentum is of order Q. Finally, integrating term by term removes the hard contribution, 
and thus we are left with the unrenormalized soft part, in which ultraviolet divergences may 
still occur (which are removed by applying R). Alternative implementations can be found 
in [28-30]. 

Since the prescription for RSG contains all the soft parts, the hard portion must be given 
by (G — RSG); the hard part involves only the large momentum contributions, including 
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some ultraviolet divergences. We can write the hard part as a series of local counterterms — a 
well-known result that is fundamental to the idea of effective field theories. Indeed, large 
momenta correspond to short distances that are tiny compared with the wavelengths of 
the external particles, so the effects can be described by a local interaction. Thus, we can 
perform the extra renormalization of absorbing the hard parts into the parameters of the 
lagrangian. 

In our case, we separate the baryon propagator into the Feynman and Density parts. 
The Feynman part contains the high-momentum contributions and also contains all the 
hard poles. The Density part contains all the soft poles describing the valence nucleons, 
which occur at lower momenta. (Note that the momenta in the meson propagators are 
always spacelike.) One can write the one-loop self-energy for the scalar two-loop integral as 

£**(*) = "f |^jA|(* -,)(?(,) 

= -f / (& A ° {k ~ q)G ' M ~ T I (& A ° (k ~ q)G ' M ' (59) 

The first term contains only the hard poles and can be absorbed by the processes outlined 
above: we expand either in powers of momentum or in powers of the mean fields (Furry's 
theorem). Then interchange orders of summation and integration; what's left are series in 
the mean fields. The second term contains the soft poles, but they are regularized by the 
theta function in the Density propagator and hence this contribution is finite. As a result, 
the procedure outlined in this work is analogous to Infrared Regularization. 



B. Power Counting 

To illustrate the order of the two-loop exchange integrals, we expand them in powers 
of momenta and pick out the dominant contributions. We can expand the isoscalar meson 
propagators as 

_L_ = \- x + X* -x z + ... / 6 Q) 
x+1 

where x = (k — q)/ms < 1 for the scalar meson and x = {k — q)jmy < 1 for the vector meson. 
This is essentially an expansion of meson exchange into contact and gradient interaction 
terms. Thus, one can pick out the leading (i.e., contact) terms 

s m „ 9 2 s f kF \kfd\k\ ^ 2 



9s 



16m| 



(A + P%) (61) 



2ir A m 2 v J E*(k) J E*(q) 



9v 
8m v 



(Pb ~ 2p|) (62) 



In the case of the pion two-loop exchange integral, the pion mass is too small for Eq. (60) 
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to converge. Therefore, we instead take the chiral limit (m n — > 0) to get the dominant term 

&«^5 + OW)- (63) 

These are the terms which set the scale for the exchange contributions. 

One can acquire the scale by comparison to terms at the mean field level; for instance, 
the terms in Eqs. (61) and (62) in ratio with the leading vector mean field term are (where 
we have used the vector meson field equation) 



2ml pB 

2 ~ ~T 

9v 2 5 



(64) 



(65) 



2m 2 v 



Pb 



Next, we consider the leading terms at normal density (k F = 1.3 fm 1 ) using Wl [2]; we 

get £f] EX /p B ~ 57 MeV, £$-ex/Pb ~ -34 MeV, and £^1 EX / p B ~ 17 MeV. These terms 
are roughly third order [y = 3) in the power counting. Observe that the self-energies are 
reproduced at the 10-40% level using only the leading-order contact term or with the chiral 
limit. To incorporate the corrections to these results, it is clearly most efficient to retain the 
full meson propagators in the exchange integrals. 

Exchange interactions are two-body, so they go like 0(p 2 ) in £. This is the same depen- 
dence as the two-body mean field terms at 0{y — 2). Because of the lack of a spin-isospin 
sum in the exchange terms, the exchange terms are numerically smaller and contribute at 
the same level as the 0{y = 3) terms. A deeper understanding of how loops fit into the 
finite-density power-counting scheme will require a look at higher-order terms in the loop 
expansion, which is currently under investigation. 

In Fig. 5, we show the individual and total two-loop contributions to the energy per 
particle plotted as a function of the density p = ps- Since the ratio Ps/pb decreases as 
Pb increases, one sees that the qualitative behavior of the scalar and vector contributions 
agrees with the leading-order estimates in Eq. (61) and (62). This is not surprising, since 
the two-loop exchange integrals are dominated by the contact terms, which have a simple 
density dependence. The total two-loop contribution is almost linearly dependent on the 
density, which occurs "by construction" with the appropriately refit coupling parameters. 
While this implies that the two-loop terms are relatively short-ranged, they are nonetheless 
nonlocal and hence they cannot be absorbed directly into the parameterization via Infrared 
Regularization. 



IV. SUMMARY 

In this work, we studied two-loop corrections to symmetric nuclear matter in a covari- 
ant effective field theory. The loop expansion gives a straightforward way to separate the 
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FIG. 5: The separate and total two-loop contributions to the energy per particle versus the density. 



short-distance physics from the long-distance physics. The former can be absorbed into coun- 
terterms already present in the effective lagrangian, and they are either removed by field 
redefinitions or fitted to empirical data. The remaining long-range exchange integrals are 
nonlocal and must be computed explicitly. They produce modest corrections to the nuclear 
binding curve and can be compensated by a small adjustment of the coupling parameters. 

Since exchange integrals in effective hadronic field theories have been studied for more 
than 30 years [35], it is important to enumerate the new features of our calculations. First, 
the QHD model studied here has its basis in a Lorentz-covariant, chiral-invariant, hadronic 
effective field theory that is tailored to the nuclear many-body problem [2] and that success- 
fully describes bulk and single-particle nuclear properties at the one- loop level [2-7, 11-14]. 
Second, using standard procedures of EFT, the loop expansion at finite density provides a 
systematic, well-defined treatment of the short- and long-range contributions to the inte- 
grals that can be extended to higher orders in loops [23, 31]. Third, when interpreted in 
the context of density functional theory, the exchange contributions to the energy introduce 
nonanalytic density dependence that is qualitatively different from that appearing in the 
mean-field theory. This should allow for an improved approximation to and parametriza- 
tion of the exact energy functional. In addition, as a criterion for discussing the size of the 
exchange contributions, one can readjust the coupling parameters to reproduce the nuclear 
matter equilibrium point and see if they remain of natural size; this was indeed the case 
with the two couplings adjusted here. The naturalness of the remaining parameters in the 
underlying EFT will be studied in future work. Finally, the size of the two-loop integrals 
(roughly third order in the mean field power counting) was determined. Fuller consideration 
of how the loop expansion fits into the power counting is left for future investigation. 
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